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Abstract 



Cold dense quark matter is in a crystalline color superconducting phase wher- 
ever pairing occurs between species of quarks with chemical potentials whose 
difference 6^ lies within an appropriate window. If the interaction between 
quarks is modeled as point-like, this window is rather narrow. We show that 
when the interaction between quarks is modeled as single-gluon exchange, 
the window widens by about a factor of ten at accessible densities and by 
much larger factors at higher density. This striking enhancement reflects the 
increasingly (1 -|- l)-dimensional nature of the physics at weaker and weaker 
coupling. Our results indicate that crystalline color superconductivity is a 
generic feature of the phase diagram of cold dense quark matter, occurring 
wherever one finds quark matter which is not in the color-flavor locked phase. 
If it occurs within the cores of compact stars, a crystalline color supercon- 
ducting region may provide a new locus for glitch phenomena. 



*Email address: adam@fnal.gov 
^Email address: krishna@ctp.niit.edu 
•^Email address: eugeneus@niit . edu 



I. INTRODUCTION 



At asymptotic densities, the ground state of QCD with three quarks with equal masses is 
expected to be the color-flavor locked (CFL) phase This phase features a condensate 

of Cooper pairs of quarks which includes ud, us, and ds pairs. Quarks of all colors and all 
flavors participate equally in the pairing, and all excitations with quark quantum numbers 
are gapped. 

The CFL phase persists for unequal quark masses, so long as the differences are not 
too large [§,0. In the absence of any interaction (and thus in the absence of CFL pair- 
ing) a quark mass difference pushes the Fermi momenta for different flavors apart, yielding 
different number densities for different flavors. In the CFL phase, however, the fact that 
the pairing energy is maximized when u, d, and s number densities are equal enforces this 
equality This means that if one imagines increasing the strange quark mass m^, all quark 
number densities remain equal until a flrst order phase transition, at which CFL pairing is 
disrupted, (some) quark number densities spring free under the accumulated tension, and 
a less symmetric state of quark matter is obtained 0. This state of matter, which results 
when pairing between two species of quarks persists even once their Fermi momenta differ, 
is a crystalline color superconductor. 

We can study crystalline color superconductivity more simply by focusing just on pairing 
between massless up and down quarks whose Fermi momenta we attempt to push apart by 
turning on a chemical potential difference, rather than a quark mass difference. That is, we 
introduce 

/i„ = yU - 5/i 

fj,d = fi + Sjj . (1) 

If 6fi is nonzero but less than some 6fii, the ground state is precisely that obtained for 6fi = 
jSHlOi. In this state, red and green up and down quarks pair, yielding four quasiparticles 



with superconducting gap Aq [^1|-|T4[. Furthermore, the number density of red and green up 
quarks is the same as that of red and green down quarks. As long as 6fi is not too large, this 
BCS state remains unchanged (and favored) because maintaining equal number densities, 
and thus coincident Fermi surfaces, maximizes the pairing and hence the gain in interaction 
energy. As 6fi is increased, the BCS state remains the ground state of the system only as 
long as its negative interaction energy offsets the large positive free energy cost associated 
with forcing the Fermi seas to deviate from their normal state distributions. In the weak 



coupling limit, in which Ao//i <^ 1, the BCS state persists for 5fj. < Sjii = Aq/v^ p|,[TO 

These conclusions are the same (as long as Ao//i <C 1) whether the interaction between 
quarks is modeled as a point-like four-fermion interaction or is approximated by single-gluon 
exchange . The latter analysis [[l5l -|2^J^,^ is of quantitative validity at densities which are 



so extremely high that the QCD coupling g{fi) < 1 ||23[. Applying these asymptotic results at 
accessible densities nevertheless yields predictions for the magnitude of the BCS gap Aq, and 
thus for 6fii, which are in qualitative agreement with those made using phenomenologically 
normalized models with simpler point-like interactions 0,^0. 

Above the BCS state, in a range < 6fi < 6fi2, the crystalline color superconducting 
phase occurs. We shall demonstrate in this paper that applying asymptotic results at ac- 
cessible densities yields a window {6^2 — ^^^i) that is more than a factor of ten wider than 
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in models with point-like interactions. The crystalline color superconductivity window in 
parameter space may therefore be much wider than previously thought, making this phase 
a generic feature of the phase diagram for cold dense quark matter. The reason for this 
qualitative increase in 6^2 can be traced back to the fact that gluon exchange at weaker 
and weaker coupling is more and more dominated by forward-scattering, while point-like 
interactions describe s-wave scattering. What is perhaps surprising is that even at quite 
large values of g, gluon exchange yields an order of magnitude increase in 6fi2 — Sfii. 

The crystalline color superconducting state is the analogue of a state first explored by 
Larkin and Ovchinnikov EM and Fulde and Ferrell P5[ in the context of electron supercon- 



ductivity in the presence of magnetic impurities. Translating LOFF's results to the case of 



interest, the authors of Ref. |jTO| found that for 6fi ^ 6^1 it is favorable to form a state in 
which the u and d Fermi momenta are given by fiu and Hd as in the absence of interactions, 
and are thus not equal, but pairing nevertheless occurs. Whereas in the BCS state, ob- 
tained for 6fi < 6 III, pairing occurs between quarks with equal and opposite momenta, when 
6fi ^ 6fii it is favorable to form a condensate of Cooper pairs with nonzero total momentum. 
This is favored because pairing quarks with momenta which are not equal and opposite gives 
rise to a region of phase space where each quark in a Cooper pair can be close to its own 
Fermi surface, even when the up and down Fermi momenta differ, and such pairs can be cre- 
ated at low cost in free energy.[] Condensates of this sort spontaneously break translational 
and rotational invariance, leading to gaps which vary periodically in a crystalline pattern. 
If in some shell within the quark matter core of a neutron star (or within a strange quark 
star) the quark chemical potentials are such that crystalline color superconductivity arises, 
as we now see occurs for a wide range of reasonable parameter values, rotational vortices 



may be pinned in this shell, making it a locus for glitch formation . Rough estimates of 
the pinning force suggest that it is comparable to that for a rotational vortex pinned in the 
inner crust of a conventional neutron star, and thus may yield glitches of phenomenological 
interest |jT(]| . 

As in Refs. |ll0| , p9| , we shall restrict our attention here to the simplest possible "crystal" 
structure, namely that in which the condensate varies like a plane wave: 

(^(x)V'(x)) oc Ae^^'i-'^ . (2) 

Wherever this condensate is favored over the homogeneous BCS condensate and over the 
state with no pairing at all, we expect that the true ground state of the system is a condensate 
which varies in space with some more complicated spatial dependence. The phonon associ- 
ated with the plane wave condensate (0) has recently been analyzed PD[. A full analysis of 



^LOFF condensates have also recently been considered in two other contexts. In QCD with 
/in < 0, /id > and Hu = —fJ-d, one has equal Fermi momenta for u antiquarks and d quarks, BCS 
pairing occurs, and consequently a (ud) condensate forms [^12?]. If — /i^ and fi^ differ, and if the 



difference lies in the appropriate range, a LOFF phase with a spatially varying (ud) condensate 
results p6| , p^ . Our conclusion that the LOFF window is much wider than previously thought 
applies in this context also. Suitably isospin asymmetric nuclear matter may also admit LOFF 



pairing, as discussed recently in Ref. [28 
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all the phonons in the crystalline color superconducting phase must await the determination 
of the actual crystal structure. 

The authors of Refs. p!0| , |29[] studied crystalline color superconductivity in a model in 
which two flavors of quarks with chemical potentials (|l|) interact via a point-like four-fermion 
interaction with the quantum numbers of single gluon exchange. In the LOFF state, each 
Cooper pair has total momentum 2q with |q| ^ 1.26fi |jTO|. The direction of q is chosen 
spontaneously. The LOFF phase is characterized by a gap parameter A and a diquark 
condensate, but not by an energy gap: the quasiparticle dispersion relations vary with the 
direction of the momentum, yielding gaps which range from zero up to a maximum of A. 
The condensate is dominated by those "ring-shaped" regions in momentum space in which 
a quark pair with total momentum 2q has both members of the pair within approximately 
A of their respective Fermi surfaces. 

The gap equation which determines A was derived in Ref. [0 using variational methods, 
along the lines of Refs. P3| , |HT| . It has since been rederived using now more familiar methods, 
namely a diagrammatic derivation of the one-loop Schwinger-Dyson equation with a Nambu- 
Gorkov propagator modified to describe the spatially varying condensate [^. This gap 
equation can then be used to show that crystalline color superconductivity is favored over 
no ud pairing for 6fi < 5^2- If quarks interact via a weak point-like four-fermion interaction, 
(5/^2 ~ 0.754Ao p^j25|j31|JT0[] . For a stronger four-fermion interaction, (5/i2 tends to decrease, 
closing the crystalline color superconductivity window ||1CI|| . 

Crystalline color superconductivity is favored for djii < 6fi < 6fi2- As 6fi increases, one 
finds a first order phase transition from the ordinary BCS phase to the crystalline color 
superconducting phase at 5/i = 6fii and then a second order phase transition at Sfi = Sfi2 
at which A decreases to zero. Because the condensation energy in the LOFF phase is much 
smaller than that of the BCS condensate at 5/1 = 0, the value of 6fii is almost identical to 
that at which the naive unpairing transition from the BCS state to the state with no pairing 
would occur if one ignored the possibility of a LOFF phase, namely = Aq/v^- For all 
practical purposes, therefore, the LOFF gap equation is not required in order to determine 
Sfii. The LOFF gap equation is used to determine Sfi2 and the properties of the crystalline 



color superconducting phase 10 



In this paper, we generalize the diagrammatic analysis of Ref. to analyze the crys- 
talline color superconducting phase which occurs when quarks interact by the exchange of a 
propagating gluon, as is quantitatively valid at asymptotically high densities. At weak cou- 
pling, quark-quark scattering by single-gluon exchange is dominated by forward scattering. 
In most scatterings, the angular positions of the quarks on their respective Fermi surfaces 
do not change much. As a consequence, the weaker the coupling the more the physics can 
be thought of as a sum of many (1 + l)-dimensional theories, with only rare large- angle 
scatterings able to connect one direction in momentum space with others [|l7l. Suppose for 
a moment that we were analyzing a truly (1 + l)-dimensional theory. The momentum-space 
geometry of the LOFF state in one spatial dimension is qualitatively different from that in 
three. Instead of Fermi surfaces, we would have only "Fermi points" at ±yU„ and i/i^. The 
only choice of |q| which allows pairing between u and d quarks at their respective Fermi 
points is |q| = 6fi. In (3 + 1) dimensions, in contrast, |q| > is favored because it allows 
LOFF pairing in "ring-shaped" regions of the Fermi surface, rather than just at antipodal 



points [p4|,p5|,p!0| . Also, in striking contrast to the (3 + l)-dimensional case, it has long 
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been known that in a true (1 + 1) -dimensional theory with a point-hke interaction between 
fermions, 8^2/^0 — 00 in the weak-interaction hmit [^. 

We therefore expect that in (3 -t- l)-dimensional QCD with the interaction given by 
single-gluon exchange, as /2 — > 00 and g{il) the |q| which characterizes the LOFF phase 
should become closer and closer to 6fi and 6^12/ Aq should diverge. We shall demonstrate 
both effects, and shall furthermore show that both are clearly in evidence already at the 
rather large coupling g = 3.43, corresponding to /2 = 400 MeV using the conventions of 
Refs. ||T9|j23[| . At this coupling, 6(12/^0 ^ 1.2, meaning that (5/i2 — ~ (1.2 — 1/a/2)Ao, 



which is much larger than (0.754 — 1/\/2)Aq. If we go to much higher densities, where the 
calculation is under quantitative control, we find an even more striking enhancement: when 
g = 0.79 we find 6^12/ Aq > 1000! We see that (relative to expectations based on experience 
with point-like interactions) the crystalline color superconductivity window is wider by more 
than four orders of magnitude at this weak coupling, and is about one order of magnitude 
wider at accessible densities if weak-coupling results are applied there. 

In Section II, we present the gap equation derived from a Schwinger-Dyson equation as 
in Ref. [^, but now assuming that quarks interact via single-gluon exchange. As our goal 
is the evaluation of 5/i2, the value of 6fi at which the crystaUine color superconducting gap 
vanishes, in Section III we take the A — limit of the gap equation. Analysis of the resulting 
equation yields 6^2- In Section IV we look ahead to implications of our central result that 
the crystalline color superconductivity window is wider by (an) order(s) of magnitude than 
previously thought. 



II. THE GAP EQUATION FOR CRYSTALLINE COLOR 
SUPERCONDUCTIVITY AT WEAK COUPLING 



In the crystalline color superconducting phase p^ , p5| , p!0| , p9| , the condensate contains pairs 
of u and d quarks with momenta such that the total momentum of each Cooper pair is given 
by 2q, with the direction of q chosen spontaneously. Such a condensate varies periodically 
in space with wavelength 7r/|q| as in (H). Wherever there is an instability towards (||), we 
expect that the true ground state will be a crystalline condensate which varies in space like 
a sum of several such plane waves with the same |q|. As in Refs. [Il0| , ^ , we focus here on 
finding the region of 6fi where an instability towards (^ occurs, leaving the determination 
of the favored crystal structure to future work. 

We begin by reviewing the field theoretic framework for the analysis of crystalline color 
superconductivity presented in Ref. [^. In order to describe pairing between u quarks 
with momentum p + q and d quarks with momentum — p + q, we must use a modified 
Nambu-Gorkov spinor defined as 



*(p,g) 



(3) 



Note that by q we mean the four- vector (0,q). The Cooper pairs have nonzero total mo- 
mentum, and the ground state condensate (|]) is static. The momentum dependence of (H) 
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is motivated by the fact that in the presence of a crystalhne color superconducting con- 
densate, anomalous propagation does not only mean picking up or losing two quarks from 
the condensate. It also means picking up or losing momentum 2q. The basis has been 
chosen so that the inverse fermion propagator in the crystalline color superconducting phase 
is diagonal in p-space and is given by 



S-\p,q) 



i>-i + iJ^dio 

-A{p, -q) 
A(p,g) 










A{p,q) 




(4) 



where A = 7oA"''7o and A = C'-f^e"'^^ is a matrix with color (a, (3) and Dirac indices 
suppressed. Note that the condensate is explicitly antisymmetric in flavor. 2p is the relative 
momentum of the quarks in a given pair and is different for different pairs. In the gap 
equation below, we shall integrate over po and p. As desired, the off-diagonal blocks describe 
anomalous propagation in the presence of a condensate of diquarks with momentum 2q. The 
choice of basis we have made is analogous to that introduced previously in the analysis of a 



crystalline quark-antiquark condensate 



We obtain the gap equation by solving the one-loop Schwinger-Dyson equation given by 



S-\k,q)-S^\k,q) 



-9 



^,TiS{p,q)T^DTB{k-p), 



(5) 



— D'^'^Sab is the gluon propagator, S is the full quark propagator, whose inverse 
), and 5*0 is the fermion propagator in the absence of interaction, given by S 



where Dj^^ 
is given by 

with A = 0. The vertices are defined as follows: 



7mA^/2 


V 














-(7mA^/2)^; 



(6) 



In previous analyses [0,^, a point-like interaction was introduced by replacing g'^D^" 
by g'^" times a constant. Here, instead, we analyze the interaction between quarks given by 
the exchange of a medium-modified gluon and thus use a gluon propagator given by 



+ 



p. 



p^ — G{p) p^ — F{p) 



(7) 



where ^ is the gauge parameter, G{p) and F{p) are functions of po and |p|, and the projectors 



are defined as follows: 



6. 



T 
Oi 



0, p. 



-9i^u + 



Pf^Pu 
p2 



P. 



(8) 



The functions F and G describe the effects of the medium on the gluon propagator. If we 
neglect the Meissner effect (that is, if we neglect the modification of F{p) and G{p) due to 
the gap A in the fermion propagator) then F{p) describes Thomas-Fermi screening and G{p) 
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describes Landau damping and they are given in the hard dense loop (HDL) approximation 



by H 



F{p) = m 



2 P 



iPl 



Wo^ (Wo 
rvo 



IPl 



IPl 



with Qo{x) = ^ log f^-^ 



-m 



IpI 



iPo 

JpL 



iPo 

JpL 



+ 



^Po 

IpI 



(9) 



where = g'^fi^/n^ is the Debye mass for Nf = 2. 

The method of solving the Schwinger-Dyson equation has been outlined in Ref. p9| 



One obtains coupled integral equations for the gap parameters Ai, which describes particle- 
particle and hole- hole pairing, and A2,3,4, which describe particle- antiparticle and antiparticle- 
antiparticle pairing. Because only Ai describes pairing of quarks near their respective Fermi 
surfaces, this is the gap parameter which describes the gap in the quasiparticle spectrum and 
is therefore the gap parameter of physical interest. For the same reason, the gap equation is 
dominated by the terms containing Ai. Upon dropping A2,3,4 and renaming Ai — >■ A, the 
gap equation derived in Ref. can be written as: 



A(A;o 



d^p 



Hpo) 



3 sin2 MM J {2nY (po + ^i) (po - ^2) 



X 



+ 



G 



+ 



{k-py-F{k-p) {k-py-G{k-p) (k-p) 



(10) 



where 



2(3{k,p) 2/^(^5^) 2pik,-p) ^P{-p,k) . 2P{k,k) . 2piP,P) 
Cp = cos — z — cos — z COS z COS z sm — z — sm 



Cr 



2 2 2 2 

cos/3(A;, —p) cos/3(— p, k) — cos P{k,p) cos/5(p, k) 



. 2l3{k,k) . 2(3{P,P) 
— 2 sm sm 



cosa(/e,p) I cosa(p, /c) sm \- cos a[—p, —k) sm — - — I 

/ , N / r ,\ ■ 2(^iP^^) , f ,^ • 2P{-P^^Y 

cosa{— k,—p) \ COS a[p^k)sm h cos a(— p, — fcj sm 



Cf = sin 



sm 



. 2/^(^5^) • 2Pi.P^P) ■ 2/5(^'-p) • 2/^(-P'^) 

— sm sm sm sm 



2/3(A;,p) . 2pi.P^^) 
2 

-|- cosq;(A;,p) cosa{p, k) siv? ^ h cosq;(— p, —k) sin^ 



2 

m.-p) 



2 

m,pY 



-f cosa(— /c, — p) cosQ;(p, /e) sm h cos a(— p, — k) sm 



fill 



and 



+ HIp + ql - lp - ql) + Iv (lp + ql + lp - ql - 2/^)' + 4A2 sin^ 



E2{p) = - - i (lp + ql - lp - ql) + W(lp + ql + lp - ql - + sin^ (12) 
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with 



cos a{k, p) = {k — q) ■ {k — p) , 

cos P{k,p) = {q + k) ■ {q - p) . (13) 

In the next section, we shall use this gap equation to obtain 5/i2, the upper boundary of the 
crystalline color superconductivity window. 

III. THE ZERO-GAP CURVE AND THE CALCULATION OF (5/i2 

Solving the full gap equation, Eq. ([T0|), is numerically challenging. Fortunately, our task 
is simpler. As 6fi 6fi2, the gap A ^ 0. Therefore, in order to determine 6fi2 we need 
only analyze the A ^ limit of Eq. (p^. Upon taking this limit, we shall find a "zero-gap 



curve" relating 6fi and |q|, as obtained for a point-like interaction in Refs. p5| , pT| , pI| . The 
largest value of 6fi on this curve is 6fi2, and the |q| at this point on the curve is the favored 
value of |q| in the crystalline color superconducting phase for 6fi 6fi2. 

To take the A — > limit, we first divide both sides of Eq. (|T^) by A^ko). We must be 
careful, however: simply cancelling the A(po)/A(A;o) which now occurs on the right-hand side 
would yield an integral which diverges in the ultraviolet. This divergence is in fact regulated 
by the po-dependence of A(po) itself. While the precise shape of A(po) remains unknown 
unless we solve the full gap equation (0), we assume A(po) has a similar po dependence as in 



the BCS case |I9| , |23| , where it is relatively flat for po < and then decreases rapidly to zero 
as po increases to several times fi. This rapid decrease regulates the integral. This means 
that for ko -C /2, we can approximate A{po)/A{ko) by a step function which imposes an 
ultraviolet cutoff on the po integral at Pq = A, where A is of order several times fi. We shall 
investigate the dependence of our results on the choice of A. After the factor A(po)/A(/co) 
has been replaced by a cutoff, it is safe to simply set A to zero everywhere else in the 
integrand. 

After rotating to Euclidean space we obtain 

d^p 1 



3 sin2 M J (27r)4 (po - «^i)(po + ^E2) 



X 



Cf _^ Co ^ 



{k-pY + F{k-p) {k-py + G{k-p) {k - pf 



(14) 



where the energies Ei^2 can be obtained from Eq. (|T^) by setting A equal to 0. We also 
simplify the functions F and G to 

F{p) = m\ G{p) = l^'f^y (15) 

valid for Po ^ |p| ~ and to leading order in perturbation theory. We shall work in the 
gauge with ^ = 0, as this will allow us to compare our results for 6^2 to values of Aq obtained 
in the same gauge in the analysis of the 6fi = |q| = version of Eq. (|r^) in Ref. []T^. As 



the gauge dependence in the calculation of Ag decreases as g is reduced below ~ 1 
we expect that the same is true here for 6fi2- We leave an analysis of the ^-dependence 
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of (5/i2 to the future. In evaluating the right-hand side of Eq. (0), we should set ko ~ A, 
but instead choose kg = for simplicity as we expect A(A;o) to be almost constant for 
ko < A [p!9|,p3|. We must choose k in the region of momentum space which dominates the 



pairing: for 5 fx = |q| = 0, for example, one takes |k| = /i |1^J2^. Now, we choose k on 
the ring in k-space which describes pairs of quarks both of whose momenta (k + q and 
k — q) lie on their respective noninteracting Fermi surfaces, as this is where pairing is most 



important [^,^,0. For |q| = 6fi ^ 0, as opposed to |q| > 7^ 0, the ring degenerates to 
the point in momentum space with |k| = /2 and k antiparallel to q. 

We do the po integral analytically. There are six poles in the complex pq plane: two from 
the quark propagator located at iEi at —iE2 and four from the gluon propagator located 



at ±z-y/ |k — pp + m?, —i^nm^ ± -y/TT^m^ + 64|k — p|^)/(8|k — p|). The poles from the gluon 
propagator have residues which are smaller than those from the quark propagator by factors 
of ft. We therefore keep only the poles at iEi and —iE2. Furthermore, we can drop the p^ 
pieces in the gluon propagators because at the E12 poles they are negligible. Upon doing 
the contour integral over po? we notice that we only obtain a non-zero result if both Ei and 
E2 are positive. This defines the "pairing region" V of Refs. |l3;29]. We get 



(16) 



3sin2M|^7p (27r)3 (Ei + E2) \\k - + F |k-p|2 + G 



where the integral is taken over the pairing region V with an upper cutoff on the p integral 
of A. The remaining integrals are done numerically. 

Upon doing the integrals, Eq. ([I6| ) becomes an equation relating |q| and 6fi. We explore 
the (|q|, 5/i) plane, evaluating the right-hand side of Eq. ([T6|) at each point, and in so doing 
map out the "zero-gap curves" along which Eq. ([T6|) is satisfied. Several zero-gap curves 
are shown in Fig. |1|, in order to exhibit the dependence of the curve on the coupling g and 
on the parameter A. The top panel is for g = 3.4328: if we use one-loop running with 
Aqcd = 200 MeV and Nf = 2, this corresponds to g{fi) for p, = 400 MeV. The lower panel is 
for g = 2.2528 corresponding in the same sense to /i = 10^ MeV. The three curves in each plot 
are drawn with A/ ft = 2.0, 2.5, and 3.0. The choice of A makes some difference to the scale of 
the zero-gap curves, but does not change the qualitative behavior. We are only interested in 
the qualitative behavior, since the couplings g which are appropriate at accessible densities 
are much too large for the analysis to be under quantitative control anyway. In subsequent 
calculations, we shall take A//i = 2.5 as our canonical choice. Removing the artificially 
introduced parameter A would require solving the full gap equation ([T0|). 

We only show the zero-gap curve for |q| > 6fi, because the favored value of |q| must 



lie in this region P^ , pB| , pT| , pi]| . Each zero-gap curve in Fig. ^ begins somewhere on the line 
Sn = \q\, moves up and to the right following the line 5/i = |q| very closely, and then turns 
around sharply and heads downward, eventually reaching the 5/1 = axis. This behavior 
is most clearly seen in the upper panel. The smaller the value of g, the more sharply the 
curve turns around and the more closely the downward-going curve hugs the 6fi = |q| line. 
Solutions to the full gap equation ( p!OD with A 7^ occur in the region bounded by the 
zero-gap curve, between it and the = |q| line. As g is decreased, this region becomes a 
sharper and sharper sliver, squeezed closer and closer to the = |q| line. Since in (1 + 1) 
dimensions |q| = is favored, the change in the zero-gap curves from the upper panel to 
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FIG. 1. Zero-gap curves for different choices of the cutoff A. The upper panel is for g = 3.4328, 
corresponding to Ji = 400 MeV, while the lower panel is for g = 2.2428, corresponding to 
fL = 10^ MeV. The three curves in each panel correspond to A//i = 2.0 (dotted red), 2.5 (solid 
black), and 3.0 (dashed blue). The square, triangle, and circle mark the locations of (5^2 for each 
of the three curves. 
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iqi/Ao 

FIG. 2. Zero-gap curves for different choices of fi. Each curve is normahzed to the BCS gap Aq 
for the corresponding fi. The three curves are p, = 400 MeV (sohd blue), 10^ MeV (dashed red), 
and 10^ MeV (dotted black). Again, 5//2 is marked for each curve. For comparison, the star is the 
position of for a weak point-like interaction. 



the lower vividly demonstrates how the physics embodied in the gap equation (|TOD becomes 
effectively (1 + l)-dimensional as g is reduced. 

The largest value of on the zero-gap curve is S^2- For 6fi < Sfi2, a range of values 
of |q| can be found for which a crystalline color superconducting phase with A 7^ exists. 
The |q| within this range which results in a phase with the lowest free energy is favored, and 
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(MeV) g 6fi2/fL H/Sf^2 Ap (MeV) Sfii/Ap 

400 3.4328 0.3317 1.012 107 1.24 

10^ 2.2528 0.1362 1.001 37.5 3.63 

10^ 1.4450 0.0241 1.000 11.9 20.3 

10^ 1.1464 5.23 xlO-3 1.0000 8.0 65.3 

10^ 0.9793 1.71 xlQ-'^ 1.0000 7.63 225 

10"^ 0.8689 5.84 xlO""^ 1.0000 9.16 638 

10^ 0.7890 1.63 xlO"^ 1.0000 13.0 1254 



TABLE I. Positions of 6fj,2 and corresponding |q| for seven values of g: if we use one- loop 
running with Aqcd = 200 MeV, these correspond to g{fl) for the values of p, shown. For each g, we 
give (5/i2/A ^-iid |q|/(^/i2 as determined from the zero-gap curves we calculate. To obtain 5^2/^0-, 
we take Aq from Ref. [|l9|. 



the corresponding A obtained by solving Eq. (^) at that |q| is favored. At 5^2-, however, 
A — >■ 0. For each curve in Fig. |l|, the location of 5^2 is marked. In all cases, |q| at 5^2 is 
very close to 5/i2- 

We see from Fig. |I| that 6^2/^ decreases as g decreases, but this is not the important 
comparison. Recall that the crystalline color superconducting phase occurs within a window 
Sfii < SjJL < 6^2 where S^i = Ao/a/2 at weak coupling, with Aq the BCS gap obtained at 
6fi = 0. To evaluate the width of the window, we must therefore also compare 6fi2 to Aq. We 



obtain Aq from Ref. |T^, and in Fig. ^ we plot the zero-gap curves for fi equal to 400 MeV, 
10^ MeV, and 10"^ MeV, in each case scaled by the corresponding BCS gap Aq. In contrast, 
the star marks the location of 6fi2 obtained in a theory with a weak point-like interaction: 
5/^2 = O.754A0 and |q| = 0.906Ao. 

The effect of the single-gluon exchange interaction, then, is to increase 6^12/^0 dramati- 
cally. Because this interaction favors forward scattering at weak coupling, it causes 6^2 to 
approach the line |q| = Sfj, and 6(12/ '^0 to diverge, as occurs in (1 -|- 1) dimensions where 
there are only "Fermi points". In Table |, we show the 8^2/^0 values for different values 
of /i, with all calculations done for A/fl = 2.5. Already aX g = 3.4328, corresponding to 
fi = 400 MeV, (5/X2 is within 1% of the |q| = 5/i line and 6fi2 = 1.24Ao. Although this result 
will depend somewhat on A, it means that the crystalline color superconducting window 
is about ten times wider than for a point-like interaction. Once g = 1.445, correspond- 
ing to /i = 10^ MeV, 6fi2 is closer than one part in a thousand to the |q| = 6fi line and 
6fi2 = 2O.3A0, meaning that the crystalline color superconducting window is about four 
hundred times wider than for a point-like interaction. 



IV. IMPLICATIONS AND FUTURE WORK 



We have found that 8^2/^0 diverges in QCD as the weak-coupling, high-density limit is 
taken. Applying results valid at asymptotically high densities to those of interest in compact 
stars, namely /2 ~ 400 MeV, we find that even here the crystalline color superconductivity 
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window is an order of magnitude wider than that obtained previously upon approximating 
the interaction between quarks as point-hke. 

This discovery has significant imphcations for the QCD phase diagram and may have 
significant imphcations for compact stars. At high enough baryon density the CFL phase 
in which all quarks pair to form a spatially uniform BCS condensate is favored. Suppose 
that as the density is lowered the nonzero strange quark mass induces the formation of some 
less symmetrically paired quark matter before the density is lowered so much that bary- 
onic matter is obtained. In this less symmetric quark matter, some quarks may yet form a 
BCS condensate. Those which do not, however, will have differing Fermi momenta. These 
will form a crystalline color superconducting phase if the differences between their Fermi 
momenta lie within the appropriate window. In QCD, the interaction between quarks is 
forward-scattering dominated and the crystalline color superconductivity window is conse- 
quently wide open. This phase is therefore generic, occurring almost anywhere there are 
some quarks which cannot form BCS pairs. Evaluating the critical temperature above 
which the crystalline condensate melts requires solving the nonzero temperature gap equa- 



tion obtained from Eq. (|T0|) as described in Ref. [^, but we expect that all but the very 
youngest compact stars are colder than Tc. This suggests that wherever quark matter which 
is not in the CFL phase occurs within a compact star, rotational vortices may be pinned 
resulting in the generation of glitches as the star spins down. 

Solidifying the implications of our results requires further work in several directions. 
First, we must confirm that pushing Fermi surfaces apart via quark mass differences has 
the same effect as pushing them apart via a Sfj, introduced by hand. Second, we must 
extend the analysis to the three flavor theory of interest. Third, we need to re-evaluate 6^2 
by comparing the crystalline color superconducting phase to a BCS phase in which spin-1 
pairing between quarks of the same flavor is allowed The results of Ref. suggest 
that this can be neglected, but conflrming this in the present case requires solving the full 
gap equation (|TDp. And, fourth, before evaluating the pinning force on a rotational vortex 
and making predictions for glitch phenomena, we need to understand which crystal structure 
is favored. 
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